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Abstract 

We study the onset of synchronous states in realistic chaotic neurons coupled by mutually in- 
hibitory chemical synapses. For the realistic parameters, namely the synaptic strength and the 
intrinsic current, this synapse introduces non-coherences in the neuronal dynamics, yet allowing 
for chaotic phase synchronization in a large range of parameters. As we increase the synaptic 
strength, the neurons undergo to a periodic state, and no chaotic complete synchronization is 
found. 
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I. INTRODUCTION 



Many neural networks relay on a balanced configuration of eletrical and chemical synapses 
for a normal functioning. The electrical synapse is usually associated to processes that 
require rapid responses, since the synaptic delay can be negleted. The chemical one is 
mediated by means of chemical transmitters and it is usually associated to processes that 
do not require rapid responses, since there is an intrinsic synaptic delay. 

Neural networks with electrical synapse, as well as the analogous linearly coupled oscilla- 
tors, have recently attracted much attention mainly because it provides a simple and clear 
scenario for the onset of synchronization. The basic idea behind this is that for interacting 
neurons and oscillators with electrical coupling, increasing the coupling strength leads to 
synchronous behavior. This relation is important since the more synchronous the oscillators 
are, the more information between themselves can be exchanged [1]. 

On the other hand, this relation to neurons coupled via chemical synapse is still unclear. 
In such a synapse, increasing its strength might change the neuronal dynamics, since the 
synapse itself is a dynamical system. This relation becomes even more complex if the 
chemical synapse is of the inhibitory type. In that case, while one neuron spikes the synapse 
forces the other neuron not to spike. 

For eletrical coupling it was found several types of synchronization in coupled chaotic 
oscillators. Complete synchronization [5], Generalized synchronization [h]]. There is a type 
of synchronization which appears for very small coupling strength the phase synchronization 
(PS) where the coupled chaotic oscillators have their absolute phase difference bounded but 
their amplitudes may be uncorrelated [9]. It was numerically seen in a variety of coupled 
oscillators jjj. These many types of synchronization were also found in neurons {4, 141 with 
eletrical synapses. In particular, PS was found in two electrically coupled neurons 12| and 



in small neural networks 
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The purpouse of this work is to analyze the inhibitory chemical synapses in coupled 
chaotic neurons, and its role for synchronization. We show that there is no complete chaotic 
synchronization, since as the coupling increases the neurons undergo to periodic states. 
This offers a great contrast to synapses of the eletrical type in which complete chaotic 
synchronization is commomly found. We also show that this inhibitory synapse is responsible 
for introducing phase synchronous behavior, for a wide range of parameters. This result is 
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biologicaly meaninful since the onset of phase synchronization provides a good enviroment 
for communication with chaotic systems since in phase synchronous states one can send 
information with low probability of errors [lj]. 

This paper is organized as follows. In Sec. [TTJ we introduce the realistic Hindmarsh-Rose 
neuron model and in Sec. 1111} we introduce the model for the inhibitory chemical synapses. 
In Sec. IIVt we show the likely ocurrance of phase synchronization in two neurons with 
chemical synapse and in Sec. [V] we give the conclusions. 



II. NEURON MODEL 



The neurons are described by the Hindmarsh-Rose model which consists of four coupled 
differential equations 
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ay + bx 2 — cx 3 — dz + I 
e-y - fx 2 - gw 
fi(-z + S(x + H)) 
v(—kw + r(y + /)) 



(1) 



This model has been shown to be realistic, since it reproduces the membrane potential of 
biological neurons 6j, and it is able to replace a biological neuron in a damaged biological 
network, restoring its natural functional activity it also reproduces a series of collective 
behaviors observed in a living neural network [8|. We integrate the Hindmarsh-Rose model 
using a Runge-Kutta of order 6 with adaptative step, and set the parameters to obtain a 
spiking/bursting dynamics. The parameters of the model are: a = 1.0,6 = 3.0, c = 1.0, 
d = 0.99,e = 1.01,/ = 5.0128,£ = 0.0278,# = 1.605, k = 0.9573,/ = 1.619, fx = 0.0021, 
v = 0.0009, r = 3.000, S = 3.966. 



III. SYNAPSE MODEL 



Each synaptic connection between the neurons is modeled by a nonlinear differential 
equation that mimics the release of neurotransmitters at the synaptic cleft and its absorption 
in the post synaptic cell [ll|. The current I syn injected in the postsynaptic cell is determined 
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by the dimensionless, scaled synaptic activation S(t). 



Isyn{t) — 9synS (t)[x(t) V rev ] 

dS 5 00 (Vj ra (t)) - S(t) 
T dt S - S^VUt)) ' 



(2) 



(3) 



where V rev is the synaptic potential, V in is the presynaptic voltage, x(t) represents the 
membrane potential of the postsynapict neuron, and r is the timescale governing receptor 
binding. is given by: 



We set the parameters of the synapse equations in order to present an inhibitory effect. 
That is done by using the following parameters: V t h = —0.80, V s i ope = 1.00,V rev = —1.58, 
and 5*o > 1. 

IV. PHASE SYNCHRONIZATION 

The condition for PS can be written as 



where fap are the phases calculated from a projection of the attractor onto appropriate 
subspaces. The neurons present a non-coherent dynamics due to the two time scales, i.e. 
bursting/spiking behavior. By non-coherent dynamics we mean that there is no clear center 
of rotation in which the trajectory spirals around and also it is not possible to define a 
Poincare section for which the trajectory crosses only once each time the neuron has a 
hyperpolarization. 

The chemical synapse introduces even more non-coherence. This happens because when 
one neuron is in a spiking behavior it inhibits the other neuron, which might hyperpolarize, 
but the neuron that has been inhibited still tries to spike. This competition generates 
more non-coherence in the phase space. As a consequence, it is rather unclear how one can 
calculate the phases for such dynamics. However, it is possible to overcome this problem 
by using the conditional Poincare map {2], which is a map of the attractor, construct by 
observing it for specific times at which events occur in one neuron. Using such technique, 
we can detect PS without actually having to measure the phase. 




(4) 



|0i - g^al < c, 



(5) 
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A. PS-sets 



The conditional Poincare map is a map of the flow. In particular, it consists in observing 
the trajectory of the neuron Mj at special times rj, with the index j = 1,2 indicating the 
two neurons. We define these times of events rj, by the following rule: 

• t{ represents the time at which the membrane potential in the neuron M% reaches a 
threshold for i-th times. 

• t\ represents the time at which the membrane potential in the neuron M\ reaches a 
threshold for i-th times. 

Then, we record the trajectory position of the neuron Mj at these times rj. As a result, we 
have a discret set of points called T>j. If T>j does not spreads over the attractor of Mj, but is 
rather localized, we say that the set T>j is a PS-set. It can be shown that PS-set implies PS 
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2j. This is so, because the difference between the time at which the i-th event happens in 
both oscillators is small, which means that the time difference \r\ — r|| < 5, with 5 being a 
small constant. As a consequence, the points in the conditional Poincare map are confined. 

In Fig. [1] we show two types of T>j set. In (A), the T>j is a PS-set. One can see 
that this set is localized, and so it does not spread over the attractor, then we have phase 
synchronization. The parameters are I = 3.12 and g syn = 0.78. In (B), is a situation where 
there is no phase synchronization, for / = 3.12 and g syn = 0.76. The set T>j spreads over the 
attractor. To have a global view of the possible behaviors in this system, in Fig. [2] we show 
the parameter space in the coordinates I x g syn . In this parameter space we depicted in color 
the parameters for which we have phase synchronization, and in white parameters for which 
either chaos with no synchronous behavior or periodic states is found. The color bar at the 
right of this figure indicates the relative area in percentage of the T>j set occupation in the 
attractor projection. To assure that we have chaotic phase synchronization we also compute 
the standard deviation of the event times. We introduce the quantity Tj = rj — rj -1 , so 
its relative standard deviation is given by q = {((Tj) 2 ) — (Tj) 2 }/{Tj), where (•) represents 
the average. In Fig. [31 we show q as a function of / and g syn . We have two kinds of 
transitions for phase synchronization. The first the neurons are chaotic and present a non 
synchronous behavior then increasing the synaptical coupling there is a transition to phase 
synchronization, fixing the current / = 3.12 the first transition happens around g syn ~ 0.77. 



FIG. 1: In black it is depicted the attractor and in red the set 2X-. In (A), the set T>j is localized with 
respect to the attractor, and the map T>j does not fulfill the whole attractor projection, occupying 
just partially this projection. This implies the presence of phase synchronization between the two 
neurons. The parameters are / = 3.12 and g syn = 0.78. In (B), the set T>j spreads over the 
attractor and fulfill it, this show that there is no PS, the parameters are / = 3.12 and g syn = 0.76 
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FIG. 2: The Parameter space / x g sync is depicted. The colored regions represent parameters in 
which the neurons present PS-sets, which imply phase synchronization. White regions represent 
either chaotic or periodic behavior. 

The second transition the neurons are in a periodic behavior and when we increase the 
coupling strength they undergo to chaos but phase synchronized. 

V. CONCLUSIONS 

We have shown that phase synchronization is a common behavior in neurons with in- 
hibitory chemical synapses for the Hindmarsh-Rose model. In addition, it is shown that 
there is no complete chaotic synchronization, since as the synapse strength increases the 
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FIG. 3: The Parameter space I x g sync is depicted. The colored regions represent the relative 
standard deviation of the burst time for one neuron. A positve relative standart deviation means 
that the neuron is still chaotic. The greather is the standart deviation the more chaotic the neuron 
is. White regions represent periodic behavior. 

neurons undergo to periodic states. This places the inhibitory chemical connection as a 
good candidate to explain information transmission processes regulated by means of phase 
synchronizat ion . 

The neurons present naturally non-coherent dynamics due to the two time scales provided 
by the bursting and spiking dynamics. As they are coupled by a chemical inhibitory coupling, 
they undergo to an even more non-coherent state. However, we could still detect the presence 
of phase synchronization using the conditional Poincare map. 
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